OSHTPY-HEP-T-02-001 
January 2002 



Yukawa Unification in SO (10) 



T. Blazek^, R. Dermisek* and S. Raby* 

Department of Physics, University of Southampton, Southampton, UK 
and Faculty of Mathematics and Physics, Comenius Univ., Bratislava, Slovakia 
* Department of Physics, The Ohio State University, 174 18th Ave., Columbus, 

Ohio 43210 



Abstract 

In simple SO(IO) SUSY GUTs the top, bottom and tau Yukawa couplings 
unify at the GUT scale. A naive renormalization group analysis, neglecting 
weak scale threshold corrections, leads to moderate agreement with the low 
energy data. However it is known that intrinsically large threshold corrections 
proportional to tan/5 ~ mt{Mz) /mi,{Mz) ~ 50 can nullify these t, b, t mass 
predictions. In this paper we turn the argument around. Instead of predicting 
fermion masses, we use the constraint of Yukawa unification and the observed 
values Mt, mb{mb), Mr to constrain SUSY parameter space. We find a 
narrow region survives for /i > with /i, M1/2 << rniQ, Aq ^ —1.9 mie 
and niiQ > 1200 GeV. Demanding Yukawa unification thus makes definite 
predictions for Higgs and sparticle masses. In particular we find a light higgs 
with mass m° = 114 ± 5 ± 3 GeV and a light stop with {m^jMiN ~ 450 
GeV and m^^ << m^^. In addition, we find a light chargino and a neutralino 
LSP. It is also significant that in this region of parameter space the SUSY 
contribution to the muon anomalous magnetic moment a^J^^^ < 16 x lO^^^''. 



1 Introduction 



Grand unification with SU (5), 5*0(10) or even partial unification with SU {4)cxSU (2)lX 
SU{2)ii explains the peculiar standard model charge assignments of quarks and leptons 
and also the observed family structure JI], |], Gauge coupling unification at a scale 
Mg ~ 3 X 10^^ GeV in supersymmetric grand unified theories [SUSY GUTs] fits the 
low energy data well 0, ^, 0. Moreover SO(IO) SUSY GUTs have many profound 
features 

• All fermions in one family sit in one irreducible 16 dimensional representation. 

• The two Higgs doublets, necessary in the minimal SUSY standard model [MSSM] , 
sit in one irreducible 10 dimensional representation. 

• Right-handed neutrinos, necessary for a see-saw mechanism for neutrino masses, 
are naturally included in the 16 dimensional representation. 

In addition, in the simplest version of SO(IO) the third generation Yukawa couplings 
are given by a single term in the superpotential = A 16 10 16 resulting in Yukawa 
unification At = A;, = = A^,^ = A. Hence, like gauge coupling unification, there is a 
prediction but this time for Mf = 180 ± 15 GeV with large tan/3 ~ 50 (see for example, 
Anderson et al. [j^); in good agreement with the data. Note, GUT scale threshold 
corrections to this Yukawa unification boundary condition are naturally small ( < 1% ), 
unlike the corrections to gauge coupling unification which can easily be several percent 
(see the Appendix for a discussion of perturbative GUT scale threshold corrections to 
gauge and Yukawa couplings). 

This beautiful prediction is however severely weakened by potentially large weak 
scale threshold corrections proportional to tan/9 [§, ^ 0. The complete set of one loop 
corrections is given by 

mb{Mz) = Xb{Mz) ^ cos(3 (1 + Ami + + ^^f + ^'^b^ + Amf^). (1) 

v2 

The first three terms are SUSY mass insertion corrections. The dominant contributions 

Aml^'^^tanP and Amf ^ ^ ^ tanP (2) 

Stt mt 1671^ mj 

t 

can be as large as 50%. Note in most regions of SUSY parameter space these two terms 
have opposite sign. Am^ is on the other hand small, 0(-l%). The log term results from 
finite wave function renormalization of the bottom quark; it is positive, independent of 
tan/? and the total from all sources is of order 6%. Finally Amf^, due to Higgs, W 
and Z exchange, is negligibly small, 0(.5%). There are similar corrections to and 
mt- The chargino corrections m^^ are also proportional to tan/5, but are significantly 
smaller than Am^ since typically we have rriiy. » mi. Finally the corrections to mt 



2 



are not proportional to tan (3. The complete set of corrections can be found in the papers 



by Rattazzi and Sarid and by Pierce et al. |T0| . 

In most regions of SUSY parameter space Am^ is dominant and in our conventions 
Am^ > for yU > 0. The sign of /i is constrained by experiment; in particular b —>■ s'-f and 
^TVBVF both favor > 0. The same one loop graphs with a photon or gluon insertion and 
outgoing b replaced with s contributes tob ^ sj. The chargino term typically dominates 
and has opposite sign to the SM and charged Higgs contributions, thus reducing the 
branching ratio for /i > 0. This is necessary to fit the data since the SM contribution is 
somewhat too big. fi < would on the other hand constructively add to the branching 
ratio and is problematic. In addition, the recent measurement of the anomalous magnetic 
moment of the muon a^^^ = (g-2)/2 = 43 (16) x 10~^° also favors /i > |n|.Q 

In a recent letter |T3[ we showed that Yukawa unification with /i > 0, including the 
complete one loop threshold corrections, is only consistent with the data in a narrow 
region of SUSY parameter space with fi, M1/2 ~ 100 — 500 GeV; Aq ~ —1.9 mie; 
mio ~ 1.4 mi6 and mig > 1200 GeV. The parameters mie, mio denote the soft 
SUSY breaking mass terms for squark and slepton, Higgs multiplets, respectively. Note 
the requirement of Yukawa unification thus dramatically constrains the SUSY particle 
spectrum and Higgs masses. 

In this paper we present a more detailed analysis of the SUSY particle spectrum 
and the allowed parameter range. In addition to fitting electroweak data and the top, 
bottom and r masses, we also include constraints from 6 — s 7 and Bs fi~^ 
The latter constraints increase the predicted stop mass and the mass of the CP odd 
Higgs A^, the heavy CP even Higgs and charged Higgs at the expense of a small 
increase in x^- We find a light CP even Higgs boson with mass = 114 ± 5 ± 3 GeV. 
In addition in the region where mig < 2000 GeV, we find a light chargino with mass 
~ 120 — 240 GeV, a neutralino LSP with mass m^^o ~ 75 — 160 GeV, a light stop with 
mass m^^ ~ 450 — 540 GeV << . The first and second generation squark and slepton 
masses are of order mig. It is also significant that in this region of parameter space we 
find a^.^^^ < 16 x 10"^°. Note also some recent discussions of Yukawa unification and 



It is well known that electroweak symmetry breaking with large tan/3 and mig >> 
Ml/2 requires Higgs up/down mass splitting [|T6|. We find however that the fits to third 
generation fermion masses are sensitive to the mechanism used to split the Higgs masses. 
In this paper we consider D term and "Just So" Higgs splitting (defined in the text). We 
study the sensitivity of our results to small GUT scale threshold corrections to Yukawa 
couplings. Significantly larger threshold corrections are needed for D term splitting 
versus the Just So case. 

The paper is organized as follows. In section |^ we discuss the analysis. We give the 
results for the case of Just So Higgs splitting in section ^ and D term Higgs splitting 
in section ^. The constraints of 6 — >^ s 7 and Bg — > fi^ /i~ and additional experimental 
tests are considered in section ^ For the impatient reader we present detailed results 



Note, recent theoretical reevaluations of the standard model contribution are now closer to expcri- 



nient with a^^^ = 25.6 (16) x IQ-^" O 
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from some typical points in SUSY parameter space in Table 1 (without) and Table 2 
(with) the constraints from 6^57 and Bg fi~^ fi~ included. Note, we have included 
the prediction for a^^^^ in Tables 1 and 2; however it has not been included in the 
function when fitting. Finally some theoretical questions are addressed in section 



2 Analysis 

We use a top - down approach with a global analysis The input parameters are 



defined by boundary conditions at the GUT scale. The 11 input parameters at Mq are 
given by — three gauge parameters Mg, ^^(Mg), €3; the Yukawa coupling A, and 7 soft 
SUSY breaking parameters /i, M1/2, tan/3; m^g, m^Q, Am|^ {Dx) for Just So (D 
term) case.0 These are fit in a global analysis defined in terms of physical low energy 
observables. We use two (one)loop renormalization group [RG] running for dimensionless 
(dimensionful) parameters from Mq to Mz- We require electroweak symmetry breaking 
using an improved Higgs potential, including and ml corrections in an effective 2 
Higgs doublet model below Mgtop IH]- The function includes 9 observables; 6 pre- 
cision electroweak data aEM^ G^, as{Mz), Mz, My/, Pnew and the 3 fermion masses 
Mtop, mb{mb), M^.g We fit the central values: Mz = 91.188 GeV, Mw = 80.419 



GeV, X 10^ = 1.1664 GeV'^, a^M = 137.04, = 1.7770 GeV with 0.1% numer- 
ical uncertainties; and the following with the experimental uncertainty in parentheses: 
as{Mz) = 0.1180 (0.0020), pnew x 10^ = -0.200 (1.1) [||, Mf = 174.3 (5.1) GeV, 
fnb{'mb) = 4.20 (0.20) GeV .0 We include the complete one loop threshold corrections at 
Mz to all observables. In addition we use one loop QED and three loop QCD RG run- 
ning below Mz- Finally, with regards to the calculated Higgs and sparticle masses, the 
neutral Higgs masses h, H, are pole masses calculated with the leading top, bottom, 
stop, sbottom loop contributions; while all other sparticle masses are running masses. 

We minimize using the CERN subroutine minuit. In order to present our results 
we typically keep three parameters (such as p, M1/2, mie) fixed and minimize with 
respect to the remaining eight parameters. We then plot our results as contours in the 
two parameter space. 



2.1 EWSB and Higgs mass splitting 

The first significant constraint derives from electroweak symmetry breaking [EWSB] in 
the large tan/5 regime. It has been shown that this typically requires m"^^ < rrijj . In 
fact more general solutions for EWSB exist with Higgs up/down splitting and with less 



fine-tuning (see |16] and Rattazzi and Sarid |^). The range of soft SUSY parameters 



required is consistent with solution (B) of Olechowski and Pokorski |T6 



^€3, defined in the Appendix, and Am^ , Dx parametrize GUT scale threshold corrections to gauge 
coupling unification and Higgs up/down mass splitting, respectively. 

■^Capital M is used for pole masses and lower case m for MS running masses. 

^The error for mb(mf,) ||l9|| appears to be quite conservative in view of recent claims to much smaller 
error bars lEll- 
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In our analysis we consider two particular Higgs splitting schemes, we refer to as Just 
So and D term splitting. In the first case the third generation squark and slepton soft 
masses are given by the universal mass parameter mig , and only Higgs masses are split: 
Just So Higgs splitting 



In this case we find Am|^ ~ 13 %. 

In the second case we assume D term splitting, i.e. that the D term for U{l)x is 
non-zero, where U{l)x is obtained in the decomposition of 5*0(10) SU{5) x U{l)x- 
In this second case, we have: D term splitting 



Here we find Amjj = 2 Dx/mj^ ~ 5%. 

Just So Higgs splitting does not at first sight appear to be as well motivated as 
D term splitting. In the Appendix we present two example mechanisms for Just So 
Higgs splitting. Here we present the most compelling argument. In 5*0(10), neutrinos 
necessarily have a Yukawa term coupling active neutrinos to the "sterile" neutrinos 
present in the 16. In fact for z/^ we have X,y^ Ur L Hu with X,y^ = Xt = Xb = Xt = X. In 
order to obtain a tau neutrino mass with ~ 0.06 eV (consistent with atmospheric 
neutrino oscillations), the "sterile" Ur must obtain a Majorana mass Mp^ > 10^^ GeV. 
Moreover, since neutrinos couple to (and not to Hii) with a fairly large Yukawa 
coupling (of order 0.7), they naturally distinguish the two Higgs multiplets. With A = 0.7 
and Mp^ = 10^^ GeV, we obtain a significant GUT threshold correction with Am'jj = .10, 
remarkably close to the value needed to fit the data. At the same time, we obtain a small 
threshold correction to Yukawa unification < 3%. (For more details, see the Appendix.) 

3 Results: Just So Higgs splitting 

Since the log corrections Arri^^ ~ 0(6%) are positive, they must be cancelled in order 
to obtain Amb< —2 % to fit mb. For fi > the gluino contribution is positive. The 
chargino contribution is typically opposite in sign to the gluino, since At runs to an 
infrared fixed point. At cx — M1/2 (see for example, Carena et al. f^). Hence in order to 
cancel the positive contribution of both the log and gluino contributions, a large negative 
chargino contribution is needed. This can be accomplished for —At > rrig and << m^. 
The first condition can be satisfied for Ao large and negative, which helps pull At away 
from its infrared fixed point. The second condition is also aided by large At. However 
in order to obtain a large enough splitting between and m^, large values of mig are 
needed. Note, that for universal scalar masses, the lightest stop is typically lighter than 
the sbottom. On the other hand, D term splitting with Dx > gives < mj. Recall 



^Ih^,h^) = mlo (1 T Am^) 



(3) 




= ^ 2Dx, 

)=^ie + Dx, 
= uiIq - 3Dx. 



(4) 
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Figure 1: contours for mig = 1500 GeV (Left) and niiQ = 2000 GeV (Right). 

Dx > is needed for electroweak symmetry breaking. As a result in the case of Just 
So boundary conditions excellent fits are obtained for the top, bottom and tau masses; 
while for D term splitting the best fits give mb(mb) > 4.59 GeV. In this section we give 
the results for the case of Just So Higgs splitting. The results for D term splitting are 
discussed in section ^ 

In Fig. 1, we show contours as a function of /i, M1/2 for mie = (1500) 2000 GeV. 
The shaded region in all the figures is excluded by the experimental bound on the 
chargino mass, my.+ > 103 GeV. The < 2 contour for mig = 1500 GeV is bounded 
by shaded < /i < 220 GeV, shaded < M1/2 < 380 GeV. For mie = 2000 GeV the region 
with < 1(2) is contained within the closed curve bounded by 150 (shaded) < fi < 
250 (380) GeV, 220 (shaded) < M1/2 < 450 (550) GeV. In Fig. 2, we show contours 
as a function of M1/2, mig for /i = 150 GeV; < 1 for > 2000 GeV. We see 
continues to decrease as mie increases. 

In Table 1 (Fits 1,2) we present the input parameters and output for two represen- 
tative points with universal squark and slepton masses mig = 1500 (2000) GeV with 
/i = 150 (200) GeV, M1/2 = 250 (350) GeV. We find reasonable fits (x^ < 3) only for 
> 1200 GeV. For mig < 1200 GeV, x^ increases rapidly. 

The bottom quark mass mb{mb) is given in Fig. 3 (Left) for mig = 2000 GeV. In Fig. 
3 (Right) we show that the fits improve with good fits extending to larger values of M1/2 
as niiQ increases (/i = 150 GeV is fixed). It should be clear that nih is the dominant pull 
on x^ as seen by comparing to the x^ contours of Fig. 1 (Right). 

In Fig. 4 we plot Ao/miQ at Mq as a function of /i, M1/2 for fixed mig = 2000 GeV. 
Good fits are obtained for Aq ps; —1.9 mig for all mig > 1200 GeV. Note, even though Aq 
is very large, the value of At at Mz is significantly smaller since it is driven to an infra-red 
quasi fixed point. [We come back to this point when we discuss vacuum stability issues 
later.] Finally reasonable fits require mio ~ 1.35 rriiQ (Fig. 5 Left) and Amjj ~ 0.13 
(Fig. 5 Right). 
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Figure 2: contours for /i = 150 GeV. 
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Figure 3: Contours of constant mb{mb) [GeV] for mig = 2000 GeV (Left) and for fi = 150 
GeV (Right). 
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Figure 4: Ao/niiQ contours for uiiq = 2000 GeV, with contours overlayed. 
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Figure 5: mio/mig contours (Left) and Amfj contours (Right) for mig = 2000 GeV, 
with contours overlayed. 
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Figure 6: Contours of constant mass insertion corrections to Srrib [%] from gluino loops 
(Left), chargino loops (Right) (with contours overlayed) and the total one loop cor- 
rection 6mh [%] (lower Center) for mig = 2000 GeV. 

The significant positive log correction to rrib is the main reason why Yukawa unifica- 
tion is only possible in a narrow region of SUSY parameter space. In order to compensate 
this, the chargino mass insertion contribution must be significantly larger than the gluino 
contribution. In Fig. 6 we give the gluino (Left), chargino (Right) mass insertion cor- 
rections to rrih and the total weak scale threshold correction drrii, (lower Center) for fixed 
mi6 = 2000 GeV as a function of /i, M1/2. Note in the region of < 1 the gluino 
(chargino) mass insertion corrections are large and of order 13 to 26 % (- 23 to -34 %), 
while the log correction is 5.6 to 6.6 %. These are the dominant corrections. The total 
SUSY correction to mb is -3 to - 4%. In the same region, the total SUSY correction to 
Mj is 7 - 8 and to is - 2 to - 4%. 

In summary, we have shown that good fits to b, t and r masses are only obtained 
in a narrow region of SUSY parameter space Aq ^ —1.9 mig, rriiQ ~ 1.35 rriiQ with 
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Figure 7: Contours of constant mass [GeV] with fixed mig = 2000 GeV. 

niiQ > 1200 GeV. A Just So Higgs splitting Am'jj ~ 0.13 is also required.^ 

This has interesting consequences for the Higgs and supersymmetric particle spec- 
trum. In Figs. 7 - 9 we give the H'^, H"^ masses. In Fig. 7, constant m^o con- 
tours are given as a function of M1/2 for fixed mig = 2000 GeV. We find, for < 1, 
mj<fi ~ 106 — 112 GeV. We show constant m^o contours for /x, M1/2 for fixed ttliq = 1500 
GeV in Fig. 8 (Left), for mig = 2000 GeV (Right) and as a function of M1/2, rriiQ for 
/i = 150 GeV in Fig. 8 (lower Center). We find, for < 1, K ~ 112 - 117 GeV. 
Finally, in Fig. 9, constant rriHo and contours are given as a function of fi, M1/2 for 
fixed mi6 = 2000 GeV. We find, for < 1, ~ 118 - 121 GeV and m| ~ 145 - 149 
GeV. 

In Figs. 10 - 13 we show constant mass contours for ti, t2, 61, and fi for fixed 
mi6 = 2000 GeV. We find for < 1, ~ 175 - 250 GeV, k ~ 630 - 800 GeV, 
bi ~ 500 — 650 GeV, fi ~ 250 — 500 GeV. Note, the upper bounds on squark and slepton 
masses increase as mie increases. Moreover, the first two generation squark and slepton 
masses are of order mig. 

Gaugino masses are smooth functions of fi, M1/2. The gluino mass is linear in M1/2 
and satisfies the empirical relation rrig = 2.5 Mi/2 + 25 GeV. In Fig. 14 we show constant 
mass contours for xf and for Xi, the LSP. We find, for x^ < 1? ~ 120 — 240 GeV and 
X? ~ 75-160 GeV. Finally the GUT scale parameters Mq, aciMc), 63, and A ~ .65 -.7 
and weak scale parameter tan /3 ~ 50 — 52 are weakly dependent on SUSY parameters. 

In the following we evaluate the sensitivity of our results to plausible threshold cor- 
rections to Yukawa unification at Mq. We consider two loop RG running of dimensionful 
parameters. We also artificially fix the CP odd Higgs mass m^o by applying appropriate 

^Note in our analysis we have fit as{Mz) — 0.118 ± 0.002. However if the central value for as{Mz) 
were to decrease to 0.116 we would obtain good fits for mb{mi,) in a larger region of SUSY parameter 
space and still have leaj < 5%. This is because both RG running and the gluino correction to the bottom 
mass are positive and proportional to as{Mz)- 
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Figure 8: Contours of constant mass [GeV] with fixed mig = 1500 GeV (Left); 
mi6 = 2000 GeV (Right) and with fixed ^ = 150 GeV as a function of M1/2, mig (lower 
Center) . 
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Figure 9: Contours of constant mass [GeV] (Left) and mass [GeV] (Right) with 
fixed mi6 = 2000 GeV. 




Figure 10: Contours of constant ti mass [GeV] for fixed mig = 2000 GeV (Left) and 
fixed /i = 150 GeV (Right) with constant contours overlayed. 
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Figure 11: Contours of constant t2 mass [GeV] for fixed mig = 2000 GeV (Left) and 
fixed /i = 150 GeV (Right) with constant contours overlayed. 




Figure 12: Contours of constant bi mass [GeV] for fixed mig = 2000 GeV (Left) and 
fixed /i = 150 GeV (Right) with constant contours overlayed. 
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Figure 13: Contours of constant t\ mass [GeV] for fixed mig = 2000 GeV (Left) and 
fixed yU = 150 GeV (Right) with constant contours overlayed. 
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Figure 14: Contours of constant Xi niass [GeV] (Left) and constant Xi mass [GeV] 
(Right) for fixed mig = 2000 GeV with constant contours overlayed. 
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Figure 15: contours as a function of e;,, for M1/2 = 300 GeV, /i = 150 GeV, mig = 
2 TeV. Just So scalar mass case. On the upper axis we give the equivalent value of 
My^ << Mg (the Majorana mass of the tau neutrino) which contributes to |eb| with 

ietr=o. 

penalties to the function. We then discuss the dependence of the Higgs spectrum 
and as a function of m^o. This will become important when considering the decay 
Bs /i^ yU~ in section |[ In addition in ^ we discuss constraints from the process 
b ^ s '-f. Both of these latter processes require a description of Yukawa matrices for the 
heaviest two families. 

3.1 Sensitivity to GUT scale threshold corrections 

In this section we check the sensitivity of our results to GUT scale threshold corrections 
to Yukawa unification. We define e?,, et by 

Aj = A (1 + ej) with i = b,t and A,- = A. (5) 

In Fig. 15 we give contours as a function of e^, for M1/2 = 300 GeV, /i = 
150 GeV, = 2000 GeV. We consider values of |ef,|, \et\ < 3%. It is clear that 
is not very sensitive to small Yukawa threshold corrections. 

The best motivated correction comes from integrating out a heavy tau neutrino 
with mass M^^. Neutrino oscillations consistent with atmospheric neutrino data sug- 
gest Mp^ ^ 10^^ GeV corresponding to a correction ei, = 2.6% with \et\ = 0. On the 
upper axis, in Fig. 15, we give the equivalent value of M^^ « Mq (the Majorana mass 
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of the tau neutrino) which contributes to |eb| with \et\ =0 (see the Appendix). In section 
^ we consider D term Higgs sphtting where threshold corrections are absolutely essential 
for reasonable fits. 

3.2 Two loop vs. one loop RGEs 

In the region of parameter space we consider, with mig > TeV, two loop RG running of 
soft SUSY breaking parameters may have significant consequences for sparticle masses 
as well as for electroweak symmetry breaking. We have checked however that a two loop 
RGE analysis for soft SUSY masses does not significantly affect our results. By this 
we mean that the same low energy results are obtained with small changes in the GUT 
scale parameters. In Table 1 (Fit 3) we present the comparison for mig = 2000 GeV, 
/i = 150 GeV, Ml/2 = 300 GeV fixed with all other input parameters varied to minimize 
using one and two loop RGEs. It is clear that the results are not significantly different 
from Fit 2. A small change in Aq is sufficient to guarantee positive squark masses 
squared. Of course, when one uses two loop RGEs for soft scalar masses consistency 
requires including one loop threshold corrections to these masses at the weak scale. We 
have not included the latter contributions in our calculations; thus we stick with the one 
loop RGE analysis from Mq to Mz for dimensionful parameters. 

3.3 dependence on m^o mass 

In the course of our analysis it became clear that there were two minima for x^'i with a 
low and high mass solution for the CP odd Higgs mass (see also Ref. [3). In order to 
make this behavior explicit we needed a way to choose particular values of m^o within 
the analysis. We accomplished this by adding a penalty to the function for any 
value of m^o outside a narrow range. Note, we then found minima of x^ for which this 
penalty vanished. 

In Fig. 16 we plot x^ as a function of m^o for fixed mie = 2000 GeV, // = 
150 GeV, Ml/2 = 300 GeV. The global minimum is at m^o = 110 GeV with the 
local minimum at m^o ~ 250 GeV with approximately 35% larger x^. The increased 
pull to x^ is mainly due to the p parameter. In Fig. 17 we plot the light Higgs mass vs. 
m^o. Note, at the minimum of x^, m/jo is a steeply rising function of m^o. However it 
quickly reaches a plateau with m/,o ~ 119 GeV. In Fig. 18 we see that tuho and mH± 
increase linearly with m^o. We conclude therefore that the non-SM Higgs masses are 
not constrained by Yukawa unification. Nevertheless, all the other predictions (the re- 
gion of SUSY parameter space, the light higgs mass and the sparticle spectrum) remain 
unchanged, i.e. independent of the non-SM Higgs masses. 
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Figure 16: ^ function of the CP odd Higgs mass m^. 




250 300 350 
niAO (GeV) 



500 



Figure 17: Light Higgs mass as a function of m^. 
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Figure 18: Heavy m^, light and charged Higgs masses as a function of m^. 

4 Results: D term Higgs splitting 

D term sphtting for Higgs up/down masses seems natural. We have thus performed a 
analysis with D term splitting. We find that Yukawa unification does not work in this 
case, i.e. the best obtained is > 5. It is easy to understand why D term splitting does 
not work. Dx > 0, needed for electroweak symmetry breaking, makes m| < m? already 
at Mg (see Eqn. ^). Therefore, in this case m^^ < m^^ and hence the gluino contribution 
to Arub dominates. 

It is important to check whether small threshold corrections to Yukawa unification 
can change this result. In this analysis we take \eb\, \et\ as large as 1O%0 in order to overlap 
with the parameter range considered in a recent paper ||l^ . In Fig. 19 we plot constant 

contours as a function of eb, et for mig = 2000 GeV, /i = 150 GeV, M1/2 = 300 GeV 
fixed and all other input parameters varied to minimize x^- Good fits are obtained with 
et ~ and ~ -7% or with et ~ —eb ~ 5%. These GUT scale corrections to Yukawa 
unification are significant. They are needed in this case for the RG evolution from Mq 
to Mz to drive mf^ << . This is because Yukawa couplings via the RGEs tend to 
drive squarks lighter, hence A* > A;, compensates for the unfavorable boundary condition 
with rrij < m\. Good fits are again obtained in the same narrow region of soft SUSY 
breaking parameters with Aq ~ —1.9 mig and mio ~ 1.35 mig, as in the case of Just 
So Higgs splitting and exact Yukawa unification (see Table 2). A few comments are in 
order. Neutrino threshold corrections give eb > 0, assuming the sterile neutrinos obtain 
a Majorana mass < Mg (see Appendix). In addition, although 7^ is SU{5) invariant, 

^We do not consider threshold corrections to Yukawa unification in excess of 3% to be smah. Such 
significant corrections would require additional physics explanations. 
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Figure 19: contours as a function of e^, et for M1/2 = 300 GeV, /i = 150 GeV, mie = 
2 TeV. D term splitting case. 

eb requires SU{5) breaking threshold corrections. 

In Table 1 (Fit 4) we present our best fit point with D term splitting and exact Yukawa 
unification. This best fit has > 5 and is thus unacceptable. Fit 5 on the other hand, is 
one point in parameter space with D term splitting and significant threshold corrections 
to Yukawa unification, i.e. = —0.08, et = 0.02. The results are comparable to 
Fit 2 with Just So Higgs splitting and exact Yukawa unification. This point with fixed 
rriiQ = 2000 GeV, fi = 150, M1/2 = 300 is similar to but not completely consistent 
with the results of Ref . [|I^ . Besides the fact that we always have larger values for the 



GUT scale Yukawa coupling than found in Ref. , we also require significantly smaller 
Yukawa threshold corrections. 



5 Experimental Tests & More Constraints 

The most unexpected result of this analysis is the constraint on the light Higgs mass. 
We find ~ 114±5±3 GeV where the first uncertainty comes from the range of SUSY 
parameters with ^ 1-5 and the second is an estimate of the theoretical uncertainties 
in our Higgs mass. Surely this prediction will be tested at either Run HI at the Tevatron 



or at LHC. We have used the analysis of Ref. ||T8[ which is a good approximation for 
rrit At/M^g^gY) < 0.5 with M^^^j^^^y) = (m|+miJ^/2. In Fig. 20 we plot this ratio. Note, 
due to our large value for At we typically find rrit At/M^^jj^y) ~ 0.7 which is somewhat 
outside the preferred range. We have also compared our results with FeynHiggsFast and 
find ours to be larger by about 3 GeV. We are not certain of the reason for this difference. 
Yukawa unification alone prefers light A^, and masses. However in the region 
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Figure 20: Contours of constant rrttAt/ M'^gjjgy) with fixed mig = 2000 GeV, as a function 

Of/i, Ml/2. 

of large tan (3 it has been shown that the process Bg fi^ /i~ provides a lower bound 
on rrij^o > 200 GeV (as pointed out in the recent works Dedes et al. and also Isidori 
and Retico [^) so that it is below the experimental upper bound B{Bs ^ fi~^ fi^) < 
2.6 X 10^^ (95% CL) |2^. [] This important new constraint is a consequence of the flavor 
violating quark- quark- A*^ couplings which result from the large threshold corrections to 
CKM angles obtained in the region of large tan/5 [^. This has only a minor impact 
on as discussed above. We find that increases by at most 40% for any ttiao less 
than ^ 350 GeV (Fig. 16). The light Higgs mass is rather insensitive to the value 
of m^io (Fig. 17); whereas m%,mH+ are linearly dependent on m^o (Fig. 18). We are 
thus not able to predict the other Higgs masses. Direct observation of A^, H^, 
may be difficult at the Tevatron, but should be possible at LHC. On the other hand, 
Bg fJ''^ fJ'~ is a significant constraint and may be discovered at the Tevatron. 

We find that the process 6 — 57 also provides significant new constraints on SUSY 
parameter space. In order to calculate 6 — > 57 we have included second family data 
(m,(2 GeV) = 110 ± 35 MeV, Mb - = 3.4 ± 0.2 GeV, Vcb = 0.0402 ± 0.0019) in 
the function in order to self-consistently obtain flavor violating SUSY contributions. 
We have used a parametrization of the Yukawa couplings at Mq which, though not 
completely general, fits the data well (see Appendix) .0 We find that the coefficient 
(jMSSM q£ order — Cf*'^ (see for example, Eqn. 9 in Ref. |]2^) with the chargino 
term dominating by a factor of about 5 over all other contributions. This is due to the 
light stop ti. In fact, 6 — > 57 is more sensitive to m^^ than mb{mb). This is because 
the amplitude depends on the inverse fourth power of the stop mass while chargino 
correction to the bottom mass depends only on the inverse second power. Fitting the 

''We thank K.S. Babu and C. Kolda for discussions. 

®Note, it has been shown (Ref. |Q) that 2nd - 3rd family flavor mixing can significantly affect the 
result for B{b ^ 57). 
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central value B{b — > 57) = 2.96 x 10"'' requires a heavier ti with {mi^)MiN ~ 450 
GeV; significantly larger than the range which provides the best fits to m;,. Nevertheless, 
ti is still the lightest squark with significant stop-sbottom splitting. The 61, fi masses 
also increase significantly. We now find mi,{m}j)MiN ~ 4.3. 

In Table 2 we present representative points which are consistent with both B{b ^ 57) 
and Bg ^ fi^ fi^ . Fits 1 and 2 are with Just So Higgs splitting, while Fits 3 and 4 are 
with D term splitting. Fits 1,3,4 (Table 2) correspond to the same values of /i, M1/2, rriiQ 
as Fits 2,4,5 (Table 1). Fits 1 to 3 have exact Yukawa coupling unification at Mq. Fit 
4 has €{, = —0.08, = 0.02. B{b — >• 57) is the strongest constraint in these fits. Better 
agreement between Yukawa unification and B{h ^ 57) is achieved with increasing mie. 
Also when we fit B{h S7) at +3cr, then we obtain slightly lower squark and slepton 
masses with the changes indicated in parentheses in Table 2. 

We have not reevaluated contours including the Bs — > /i^ and h ^ con- 
straints. We do not expect the contours in the /i, M1/2 plane to change significantly, 
since we can accomodate these new constraints with small changes in the parameter Aq 
and negligible changes in all others. Thus we expect that the predictions for gaugino 
masses to be unaffected by the Bg /i~ and b s'-f constraints. Hence the lightest 

neutralino is the LSP and a dark matter candidate |26||. In order to know how observable 



neutralinos and charginos may be, we encourage the analysis of some new benchmark 
points consistent with Yukawa unification. 

Finally, we recall that proton decay experiments prefer values of mig > 2000 GeV and 
niiQ » Ml/2 (see Ref. [^). This is in accord with the range of SUSY parameters found 
consistent with third generation Yukawa unification. There is however one experimental 
result which is not consistent with either Yukawa unification or proton decay and that 
is the anomalous magnetic moment of the muon. Large values of mig > 1200 GeV 
lead to very small values for a^^^ < 16 x 10" ^''.^ Hence a necessary condition for 
Yukawa unification is that forthcoming BNL data |Tl|] and/or a reanalysis of the strong 



interaction contributions to aj^*'^ will significantly decrease the discrepancy between the 
data and the standard model value of a^. 

In summary, most of the results of our analysis including only third generation 
fermions remain intact when incorporating flavor mixing. The light Higgs mass and most 
sparticle masses receive only small corrections. The lightest stop mass increases, due to 
6 ^ S7. Nevertheless there is still a significant ti — ^2 splitting and m^^ << . The 
A^, H^, masses are necessarily larger in order to be consistent with Bs /i"*^ /x" 



which suggests that this process should be observed soon; possibly at Run HI of the 
Tevatron. Finally, the central value for a^^'^ must significantly decrease. The "smok- 
ing guns" of SO(IO) Yukawa unification, presented in this paper, should be observable 
at Run HI of the Tevatron or at LHC. Also, in less than a year we should have more 



information on a^^'^ . 



Although afy^^ has not been included in the function, we have included the calculated values 



for a^^^^ for the points in Tables 1 and 2. 
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6 Discussion 



In the previous section we presented some experimental tests of Yukawa unification. Here 
we consider some open theoretical questions. 

Yukawa unification only works in a narrow region of SUSY parameter space with 
Aq ^ —1.9 mi6, mio ^ 1.35 mig and mig > 1200 GeV. The question arises, is this 
boundary condition natural in any SUSY breaking scheme? In mSUGRA, dilaton and 
anomaly mediated SUSY breaking schemes Aq ^ at the GUT scale. On the other hand, 
in gauge mediated or gaugino mediated SUSY breaking schemes, = at Mq- The 
latter are thus disfavored by Yukawa unification. In addition, anomaly mediated SUSY 
breaking has other problems. Slepton masses squared are negative unless other physics is 
added. More importantly, however, since the gluino and chargino masses have opposite 
sign, it is difficult to simultaneously fit 6 —> 57 and a^^^ . Finally although 7^ at 
Mg in mSUGRA and dilaton SUSY breaking schemes, this still does not explain why 
Aq ^ —1.9 TTiig. However, the other relations for Higgs masses, i.e. mio ~ 1.35 mie and 
Am'jj = .13, may be obtained via RG running above Mq or via threshold corrections. 

It is an interesting, but not too surprising, result that the region of SUSY parameter 
space preferred by Yukawa unification is very similar to the region of SUSY parameter 
space needed to obtain heavy 1st and 2nd generation squarks and sleptons with third 
generation squarks and sleptons lighter than O(TeV) First of all, the S'O(IO) 

boundary conditions with m^g a d l e) ~ "^i6 '^\hu hJ ~ ""^lo obtained as a 
result of demanding the inverted scalar hierarchy in Ref. [^, whereas for us they are 
input. In addition, we need a light stop with large stop-sbottom splitting forcing us to 
the same region of parameter space with large negative Aq ~ —2 mig and large mig with 
mio ~ V^miQ. 

It would be interesting to see how sensitive our results may be to alternative elec- 
troweak symmetry breaking approximations. In this paper we have used the effective 2 
Higgs doublet analysis of |TE]. This approximation may be particularly well suited to 
the light Higgs spectrum we obtain in our analysis. The alternate scheme, in which the 



Higgs tadpoles are evaluated in the MSSM at a scale of order Mstop [|T9' however 
more frequently used in the literature. 

We have a neutralino LSP and it is important to know if it is consistent with cos- 
mology and possibly a good dark matter candidate. We note that Yukawa unification 
places us in a region of SUSY parameter space which is markedly different than has been 
studied in the literature. A preliminary investigation suggests that there are no major 



problems [^. Further study in this region of parameter space is now highly motivated 
by our results. 

Finally, let us consider the issue of vacuum stability. Since we have large At, we 
may find that the vacuum stability condition \At\ < m^^ + m^^ |^ in the stop-Higgs 
sector is violated. Indeed we find only a narrow region in SUSY parameter space with 

< 1 where this constraint is satisfied. However the small change in Aq, necessary to 

""^^Note, this region of parameter space is desirable for suppressing flavor and CP violating SUSY 
loops. 
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fit B{b — s> 57), is also sufficient to satisfy tliis stability constraint. 
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A Appendix 

In this Appendix we consider the possible GUT scale threshold corrections to gauge and 
Yukawa unification and to soft SUSY breaking scalar masses. 

A.l Gauge coupling unification 

At tree level the three gauge couplings unify at the GUT scale. At one loop all three 
gauge couplings receive corrections depending logarithmically on an arbitrary scale /x 
and the masses of the particles integrated out of the theory. We may choose fi = Mq 
such that = Q!2(At) = etc- Then the one loop threshold correction corresponds to a 

shift given by = [a^lMc) — aG)/aG- ^3 obtains contributions from all massive states 
with SO(IO) quantum numbers. We obtain 



€3 =/(Cl,---,Cm) + ^log 

OTT 



detM. 



McdetM^ 



+ ■■■ (6) 



where the ffist term represents the contributions from Wsym breaking- It is only a func- 
tion of U(l) and R invariant products of powers of vevs {Q}- It is typically large 
O(±10%). The second term comes from the Higgs sector where the color triplet and 
doublet mass matrices Mt, only include those states, from 5s and 5s of SU(5) con- 
tained in Wgymbreaking and Wniggs, which mix with the Higgs sector. For further details. 



see Ref. 100 



A. 2 Yukawa coupling unification 

The third generation Yukawa couplings are derived from the minimal interaction W = 
A I63 IOh I63. At tree level we have Xt = Xt = = X^^ = X. At one loop this 



relation is corrected However in this case there are just three sources of corrections: 
gauge exchange for IOh and I63; Yukawa exchange with color triplet Higgs fields and 
with heavy right-handed neutrinos in the loop. When one considers a theory of fermion 
masses for three families then additional Yukawa couplings mixing I63 with other heavy 
SO(IO) states are possible. However even if these new Yukawa couplings are of order one, 
their contribution to threshold corrections are typically less than 1%. More importantly 
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however there are corrections which come at tree level in the effective low energy field 
theory. In a three family model 3x3 Yukawa matrices are needed in order to obtain 
both fermion masses and CKM mixing angles. Upon diagonalizing the Yukawa matrices, 
one effectively obtains tree level corrections to the simple SO(IO) relation. 



A. 2.1 One loop corrections neglecting 

Consider first the one loop corrections neglecting Vr- We find 



where Mio, M5 (M^) is the mass of the S'O(IO) gauge fields contained in S'O(IO) / SU (5), 
SU (h) / {SU {?>) X SU{2) X U{1)) (the mass of the color triplet Higgs fields contained in 
IOh). For Mr,^ Mt Mq and < IOMq we obtain eb ~ -0.7%, < 1.1%. 



A. 2. 2 One loop corrections due to ly^ alone 

The superpotential for the neutrino sector is given by 

W = XH^LrUr + ^Mp^ Ur (9) 

where Mp^ is the effective Majorana mass for the right handed tau neutrino. Neutrino os- 
cillations consistent with atmospheric neutrino data suggest Mp^ ^ mf / V3.5 x 10~^ eV^ ~ 
10^3 GeV or MpjMc ^ lO'^. 

Integrating out the right handed tau neutrino leads to equal finite wave function 
renormalization of if^, L^. Hence 

or eb = 2.6% for A = 0.7, Mg = 3 x 10^^ GeV, Mp^ = 10" GeV. This is a considerable 
correction which actually goes in the wrong direction for fitting the third generation 
masses. We discuss the sensitivity of our results to such a correction in the text. 



A. 2.3 Tree level corrections due to more realistic Yukawa matrices 



A significant threshold correction to Yukawa unification can come when Yukawa matrices 
for three families are considered. As an example, in Table 2 the process 6^57 was 
calculated using the following ansatze for two family Yukawa matrices, since it has been 



shown that they provide a good fit for fermion masses and mixing angles |32 
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^<i= \ ^' N ^ (11) 

Y.- ( 

The universal Yukawa coupling A plus the three new complex parameters e, e' and a 
and one real parameter r were varied to minimize the function with five additional 
observables = 105.66 ± 1.1 MeV, M^- M^ = 3.4 ± 0.2 GeV, m,(2 GeV) = 110 ± 35 
MeV, Vch = 0.0402 ±0.0019 and B{b s-f) x 10^ = 0.296 ±0.035. Clearly there are more 
parameters than observables, but we are not attempting in this analysis to make any 
new predictions for fermion masses. We just want to be able to calculate the branching 
ratio for 6 — > s 7 self-consistently. With Fit 1 in Table 2 as a guide, we find (upon 
diagonalizing the 2x2 Yukawa matrices at the GUT scale) ~ ~ —0.08. Several 
points should be made here. The first is that this tree level correction is model dependent 
and much larger than any one loop correction. Secondly, 7^ is a consequence of the 
Georgi-Jarlskog like mass matrices distinguishing quarks and leptons. It is needed to 
obtain a reasonable fit to ms{2 GeV) and M^. Finally, since Vet ~ 0.04 is small, we 
find et ~ efe. This latter result is in the wrong direction for obtaining good fits to third 
generation masses with D term splitting as is evident in Fig. 19 and Fit 4 (Table 2). 



A. 3 Higgs mass splitting 

D term splitting of the Hu, Hd masses is quite natural in 5*0(10) SUSY GUTs. It can 
be generated in the process of breaking 5*0(10) SU{5) by a mismatch in the vacuum 
expectation values of the 16 and 16 which are needed to break 50(10) — >■ SU{5) and 
reduce the rank of the group. Once this D term is generated it then gives mass to scalars 
proportional to their U{l)x charge. 

The Just So case does not at first sight appear to be similarly well motivated. In 
this Appendix we attempt to rectify this apparent difficulty. It is quite clear that in 
any SUSY model the Higgs bosons are very special. R parity is used to distinguish 
Higgs from squarks and sleptons. In addition, a supersymmetric mass term /i with value 
of order the weak scale is needed for the Higgs bosons. Since /i is naturally of order 
Mq, one needs some symmetry argument why it is suppressed. Of course, if the Higgs 
are special, then perhaps this will help us understand how to obtain splitting of the 
Higgs up/down while maintaining universal squark and slepton masses. GUT threshold 
corrections to soft SUSY breaking scalar masses have been considered previously. In 
Murayama et al. it was shown that the necessary condition, mio > miQ, can naturally 
be obtained in 50(10) with RG running from Mpi to Mq- In the paper by Polonsky 
and Pomarol [|1^] the splitting of the soft masses of SU (5) multiplets within irreducible 
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so (10) representations was considered. In the following we consider two novel sources 
for Higgs up/down splitting in the context of SO (10). 

A. 3.1 Ur contribution to Higgs splitting 

In the MSSM superpotential below Mq we have the contribution which distinguishes 
Hu and Hd (see Eqn. This leads to a significant threshold correction 

Am^^ ~ (2m^g + + Al) log(^) + non log terms (12) 

Using the values A = 0.7, Mp^ = lO^^ GeV and Mg = 3 x 10^^ GeV and the typical 
boundary conditions Aq ^ 2m^Q ^ 4m^g, we obtain Am|^ = ^ Am^^/m^Q = .10. Note 
this is remarkably close to the value needed for Just So Higgs splitting (see Fig. 5 (Right) 
and Fits 1 - 3 (Table 1) and Fits 1,2 (Table 2)). 

A. 3. 2 Another possible source for Higgs splitting 

Consider also the possible superpotential for the Higgs sector 

W = XaW 45 10' + X [lO'y + tfjtp 10' + M^Pij + Tr{45'^Ml^) (13) 

The first two terms are necessary for Higgs doublet/triplet splitting with assumed vevs 
for the adjoint field (45) ~ (S - L) Mg and {X) « Mq. The fields -0, are in 16,16 
dimensional representations. Their vevs break 5*0(10) to SU{b) and split the masses 
of 5io', 5io'. With this splitting and also an assumed SU{b) invariant splitting in the 
supersymmetric masses for the 45 we find significant threshold corrections to Am'jj . 
Schematically, we find 

Am^ cc ^ logiMl.JMlJ szn{e + e')sin{e - e'){r} + ■■■ (14) 

The factor r represents the ratio of soft scalar masses m\^/m\Q, m^Qz/m^Q, Ag/m^g 
and finally the dots represent the contribution of color triplet states to the loops. The 
correct sign for Arn^j can always be obtained. If A^/47r ~ 0(1) or r >> 1, then we can 
easily obtain Amjj ~ 20 — 30%. In the latter case, top, bottom unification will have 
considerably smaller threshold corrections. 
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Table 1: Five representative points of the fits. The first three are with Just So Higgs splitting 
and the last two are with D term splitting. All fits assume exact Yukawa unification at the 
GUT scale, except for 5 which has £{, = —0.08, et = 0.02. All fits are with one loop RG running 
from Mg to Mz for dimensionful parameters, except for 3 which uses two loop running. All 
entries are in units of GeV to the appropriate power. 
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51.2 


50.5 


50.6 


54.3 


51.1 


Ao/rriie 




-1.83 


-1.87 


-1.83 


-0.37 


-1.87 


observables 


Exp {(j) 












Mz 


91.188 (0.091) 


91.13 


91.14 


91.14 


91.15 


91.15 


Mw 


80.419 (0.080) 


80.45 


80.45 


80.44 


80.44 


80.44 




1.1664 (0.0012) 


1.166 


1.166 


1.166 


1.166 


1.166 




137.04 (0.14) 


137.0 


137.0 


137.0 


137.0 


137.0 


as{Mz) 


0.118 (0.002) 


0.1175 


0.1176 


0.1175 


0.1161 


0.1179 


Pnew X 10 


-0.200 (1.10) 


0.696 


0.460 


0.437 


0.035 


0.265 


Mt 


174.3 (5.1) 


175.5 


174.6 


174.4 


177.9 


174.1 




4.20 (0.20) 


4.28 


4.27 


4.28 


4.59 


4.22 


Mr 


1.7770 (0.0018) 


1.777 


1.777 


1.777 


1.777 


1.777 


TOTAL 




1.50 


0.87 


0.91 


5.42 


0.45 






116 


116 


117 


115 


116 






120 


121 


121 


117 


120 






110 


110 


111 


110 


110 


H+ 




148 


148 


149 


146 


148 


xl 




86 


130 


130 


86 


99 


xl 




135 


190 


189 


126 


152 


xt 




123 


178 


177 


105 


131 


9 




661 


913 


898 


902 


787 


ii 




135 


222 


235 


1020 


141 


bi 




433 


588 


589 


879 


342 


fi 




288 


420 


542 


1173 


693 


^SUSY X IQlU 


25.6 (16) 


9.7 


5.5 


5.5 


6.1 


6.4 



Table 2: Four representative points of the fits. The first two are with Just So Higgs sphtting 
and the last two are with D term splitting. All fits assume exact Yukawa unification at the 
GUT scale, except for 4 which has = —0.08, et = 0.02. All fits are with one loop RG 
running from Mq to Mz for dimensionful parameters. For these fits the branching ratio for 
6 — > S7 is included in y^. In addition, the CP odd Higgs mass m\ is constrained to be 200 
GeV. These points are thus consistent with both B(b ^ 57) and Bs /x"*" iJ,~. In Fits 1, 2 we 
show the change in the ii, bi, fi masses if we fit 6 —> s 7 at the central value + 3a. 



Data points 




1 


2 


3 


4 


Input parameters 
















24.72 


24.78 


24.75 


24.62 


Mg X 10-16 




3.00 


3.12 


3.09 


3.20 


es 




-0.040 


-0.041 


-0.045 


-0.043 


A 




0.63 


0.61 


0.80 


0.60 


mi6 




2000 


3200 


2000 


2000 


mio/mi6 




1.32 


1.30 


1.19 


1.30 




0.13 


0.13 


0.07 


0.05 


Ml/2 




350 


350 


350 


300 


fi 




200 


150 


115 


150 


tan P 




52.5 


50.6 


55.0 


50.7 






-1.71 


-1.83 


-0.05 


-1.75 


observables 


Exp (a) 










Mz 


91.188 (0.091) 


91.18 


91.18 


91.16 


91.18 


Mw 


80.419 (0.080) 


80.42 


80.42 


80.43 


80.42 


X 10^ 


1.1664 (0.0012) 


1.166 


1.166 


1.166 


1.166 




137.04 (0.14) 


137.0 


137.0 


137.0 


137.0 


as{Mz) 


0.118 (0.002) 


0.1172 


0.1173 


0.1162 


0.1167 


Pnew X 10 


-0.200 (1.10) 


0.228 


0.321 


0.221 


0.279 


Mt 


174.3 (5.1) 


173.8 


172.1 


178.8 


172.6 


mb{mb) 


4.20 (0.20) 


4.46 


4.42 


4.56 


4.50 




1.7770 (0.0018) 


1.777 


1.777 


1.777 


1.777 


TOTAL 




2.05 


1.72 


5.03 


3.38 






118 


119 


115 


118 






217 


217 


218 


216 






200 


200 


200 


200 


H+ 




229 


229 


228 


228 


X? 




130 


110 


86 


99 


xl 




190 


160 


126 


152 


xt 




178 


136 


105 


131 


9 




909 


904 


902 


781 


ii 




509 (-30) 


511 (-27) 


1067 


443 


bi 




749 (-42) 


903 (-14) 


900 


550 






459 (+47) 


1001 (-25) 


1173 


854 


B{h 57) X 10^ 


0.296 (0.035) 


0.297 


0.297 


0.297 


0.297 




25.6 (16) 


5.8 


2.2 


6.4 


6.4 



